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ABSTRACT 

In the standard scenario for galaxy evolution young star-forming galaxies transform into red bulge-dominated spheroids, where star 
formation has been quenched. To explain such a transformation, a strong negative feedback generated by accretion onto a central 
super-massive black hole is often invoked. The depletion of gas resulting from quasar-driven outflows should eventually stop star- 
formation across the host galaxy and lead the black hole to "suicide" for starvation. Direct observational evidence for a major quasar 
feedback onto the host galaxy is still missing, since outflows previously observed in quasars are generally associated with the ionized 
component of the gas, which only accounts for a minor fraction of the total gas content, and typically occurring in the central regions. 
We used the IRAM PdB Interferometer to observe the CO(l-O) transition in Mrk 231, the closest quasar known. Thanks to the wide 
band we detect broad wings of the CO line, with velocities up to 750 km/s and spatially resolved on the kpc scale. Such broad CO 
wings trace a giant molecular outflow of about 700 Mg/year, far larger than the ongoing star-formation rate (~200 Mg/year) observed 
in the host galaxy. This wind will totally expel the cold gas reservoir in Mrk 231 in about 10' yrs, therefore halting the star-formation 
activity on the same timescale. The inferred kinetic energy in the molecular outflow is ~ 1.2 x lO"** erg/s, corresponding to a few 
percent of the AGN bolometric luminosity, which is very close to the fraction expected by models ascribing quasar feedback to 
highly supersonic shocks generated by radiatively accelerated nuclear winds. Instead, the contribution by the SNe associated with the 
starburst fall short by orders of magnitude to account for the kinetic energy observed in the outflow. The direct observational evidence 
for quasar feedback reported here provides solid support to the scenarios ascribing the observed properties of local massive galaxies 
to quasar-induced large scale winds. 
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, 1. Introduction 

In the standard scenario for galaxy evolution emerging from both 
observations and models young star-forming galaxies transform 
into red bulge-dominated spheroids, where star-formation has 
been quenched (Bell et al. 2004). Three main drivers are pro- 
posed to explain such evolution: conditions at the time of galaxy 
formation, galaxy interactions and mergers, and AGN activity. 
The ubiquitous discovery of super-massive back holes (SMBHs) 
' at the centre of local bulges and the correlations between their 
masses and bulge properties like mass, luminosity and veloc- 
ity dispersion (Ferrarese & Ford 2005 and references therein), 
[ suggest tight links between the evolution of AGN and their host 
. galaxies. Indeed, massive galaxies not only live in reach environ- 
ments (see e.g. Bolzonella et al. 2010, Drory et al. 2010) but are 
also formed in biased regions of the density field. Furthermore, 
galaxy mergers are more common in rich environments, and 
are thought to destabilize cold gas, therefore enhancing star- 
formation and also funneling gas into the nuclear region. Such 
gas eventually accretes onto the nuclear SMBH, triggering an 
Active Galactic Nucleus (AGN). The gas and dust can intercept 
the line of sight to the nucleus, and therefore a natural expec- 
tation is that the early, powerful AGN phase is also highly ob- 
scured. Once SMBHs reach masses > 10^ Mq the AGN power 
is high enough to efficiently heat the gas and expel it from the 
galaxy through powerful winds (Silk & Rees 1998, Fabian et al. 
1999). The depletion of cold gas quenches star-formation and 



makes the QSO nucleus "suicide", forcing galaxies toward the 
BH-spheroid mass relation observed locally. Such AGN nega- 
tive feedback on the star formation in the host galaxy is thought 
to help solving two long-standing problems of galaxy evolution 
scenarios: the observed small number of massive galaxies rel- 
ative to the prediction of models (without AGN feedback) and 
their red colors, indicative of old stellar populations (Granato et 
al. 2004, Di Matteo et al. 2005, Menci et al. 2006, Menci et al. 
2008, Bower et al. 2006). 

While observational evidence for feedback on the Intra 
Cluster Medium has been observed in radio loud AGN (Fabian et 
al. 2003; McNamanara & Nulsen 2007), direct observational ev- 
idence for AGN feedback onto the gas in the host galaxy (out of 
which stars form) is still missing. In QSOs and nearby Seyfert 
nuclei prominent outflows are observed, but generally only in 
the ionized gas component (e.g. Crenshaw et al. 2003; Turnshek 
1984; Pounds et al. 2003), which accounts only for a minor frac- 
tion of the gas mass in the host galaxy, and are generally con- 
fined in the nuclear region on pc scales or in the photoionization 
cones. High-velocity energetic winds, whose kinetic energy is 
comparable to the bolometric energy of the quasar, are observed 
in the galactic nuclei (Reeves et al. 2009, Moe et al. 2009, Dunn 
et al. 2010, Bautista et al. 2010). Morganti et al. (2010) reported 
evidence for AGN induced massive and fast outflows of neutral 
H in powerful radio galaxies, possibly driven by the AGN jets. 

The bulk of the gas in QSO hosts, i.e. the molecular phase, 
appears little affected by the presence of the nuclear AGN. 
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Indeed, most studies of the molecular gas in the host galax- 
ies of QSOs and Seyfert galaxies have found narrow CO lines 
(width of a few 100 km/s), generally tracing regular rotation 
patterns, with no clear evidence for prominent molecular out- 
flows (Downes & Solomon 1998, Wilson et al. 2008, Scovifle 
et al. 2003), even in the most powerful quasars at high redshift 
(Solomon & Vanden Bout 2005, Omont 2007). Yet, most of the 
past CO observations were obtained with relatively narrow band- 
widths, which may have prevented the detection of broad wings 
of the CO lines possibly associated with molecular outflows. 
Even worse, many CO surveys were performed with single dish, 
where broad CO wings may have been confused with baseline 
instabilities and subtracted away along with the continuum. 

In this paper we present new CO(l-O) observations of Mrk 
231 obtained with the IRAM Plateau de Bure Interferometer 
(PdBI). Mrk 231 is the nearest example of a quasar object and it 
is the most luminous Ultra-Luminous Infrared Galaxy (ULIRG) 
in the local Universe (Sanders et al. 1988) with an infrared lumi- 
nosity of 3.6x10^^ Lq (assuming a distance of 186 Mpc). A sig- 
nificant fraction (~70%) of its bolometric luminosity is ascribed 
to starburst activity (Lonsdale et al. 2003). Radio, millimeter and 
near-IR observations suggest that the starbursting disk is nearly 
face-on (Downes & Solomon 1998, Cai-ifli et al. 1998, Taylor et 
al. 1999). In particular, past CO(l-O) and (2-1) IRAM PdBI ob- 
servations of Mrk 23 1 show evidence for a regular rotation pat- 
tern and a relatively narrow profile (Downes & Solomon 1998), 
and a molecular disk (Carilli et al. 1998). The existence of a 
quasar-like nucleus in Mrk 23 1 has been unambiguously demon- 
strated by observations canied out at different wavelengths, that 
have revealed the presence of a central compact radio core plus 
pc-scale jets (Ulvestad et al. 1999), broad optical emission lines 
(Lipari et al. 2009) in the nuclear spectrum, and a hard X-ray 
(2-10 keV) luminosity of lO'*'* erg s"' (Braito et al. 2004). In ad- 
dition, both optical and X-ray data have revealed that our line 
of sight to the active nucleus is heavily obscured, with a mea- 
sured hydrogen column as high as N// = 2 x 10^"* cm"^ (Braito 
et al. 2004). Mrk 231 displays clear evidence of powerful ion- 
ized outflows by the multiple broad absorption lines (BAL) sys- 
tems seen all over its UV and optical spectrum. In particular, 
Mrk 231 is classified as a low-ionization BAL QSOs, a very 
rare subclass (~10% of the entire population) of BAL QSOs 
characterized by weak [OIII] emission, in which the covering 
factor of the absorbing outflowing material may be near unity 
(Boroson & Meyers 1992). Furthermore, giant bubbles and ex- 
panding shells on kpc-scale are visible in deep HST imaging 
(Lipari et al. 2009). Recent observations with the Herschel Space 
Observatory have revealed a molecular component of the out- 
flow, as traced by H2O and OH molecular absorption features 
(Fischer et al. 2010), but the lack of spatial information has pre- 
vented an assessment of the outflow rate. 

2. Data 

We have exploited the wide bandwidth offered by the PdBI to ob- 
serve the CO(l-O) transition in Mrk 231. The observations were 
candied out between June and November 2009 with the PdBI, us- 
ing five of the 15 m antennas of the array. We observed the CO(l- 
0) rotational transition, whose rest frequency of 1 15.271 GHz is 
redshifted to 110.607 GHz (z = 0.04217), by using using both 
the C and D antenna configurations. The spectral con^elator was 
configured to cover a band width of about 1 GHz in dual polar- 
ization. The on-source integration time was ~ 20 hours. The data 
were reduced, calibrated channel by channel and analyzed by us- 
ing the CLIC and MAPPING packages of the GILDAS software. 
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Fig. 1. Continuum subtracted spectrum of the CO( 1-0) transition 
in Mrk 23 1 . The spectrum was extracted from a region twice the 
beam size (Full Width Half Maximum, FWHM), and the level 
of the underlying continuum emission was estimated from the 
region with v > 800 km s ' and v < -800 km s Left panel: 
full flux scale. Right panel: expanded flux scale to highlight the 
broad wings. The line profile has been fitted with a Gaussian 
nan^ow core (black dotted line) and a Gaussian broad component 
( long-dashed line). The FWHM of the core component is 180 
km s ' while the FWHM of the broad component is 870 km s ', 
and reaching a Full Width Zero Intensity (FWZI) of 1500 km 
s-i. 

The absolute flux was calibrated on MWC349 (S(3mm) = 1.27 
Jy) and 1150+497 (S(3mm) = 0.50 Jy). The absolute flux cal- 
ibration error is of the order +10%. All maps and spectra are 
continuum subtracted, the continuum emission being estimated 
in the spectral regions with velocity v > 800 km s ' and v < -800 
km s"'. 

3. Results 

Fig. 1 shows the spectrum of the CO(l-O) emission line, domi- 
nated by a narrow component (Full Width Half Maximum ~ 200 
km s"'), already detected in previous observations (Downes & 
Solomon 1998, Bryant & Scoville 1997). However, our new data 
reveal, for the first time, the presence of broad wings extend- 
ing to about +750 km s"', which have been missed, or possibly 
confused with the underlying continuum, in previous narrower 
bandwidth observations. Both the blue and red CO(l-O) wings 
appear spatially resolved, as illustrated by their maps (Fig. 2). 
The peak of the blue wing emission is not offset relative to the 
peak of the red wing, indicating that these wings are not due to 
the rotation of an inclined disk, leaving outflowing molecular gas 
as the only viable explanation. A Gaussian fit of the spatial pro- 
file of the blue and red wings (by also accounting for the beam 
broadening), indicates that the out-flowing medium extends over 
a region of about 0.6 kpc (0.7") in radius. To quantify the sig- 
nificance of the spatial extension of the high velocity, outflow- 
ing gas, we fitted the visibilities in the Mv-plane. We averaged 
the visibilities of the red and blue wings in the velocity ranges 
500 + 800 km s"' and -500 + -700 km sr\ and we have fitted 
a point source, a circular gaussian and an inclined disk model. 
The results of the Mv-plane fitting are shown in Fig. 3 and sum- 
marized in Table 1. The upper panels of Fig. 3 show the maps 
of the residuals after fitting a point source model. The residuals 
of the red wing are 5cr above the average rms of the map and 
those of the blue wing 3cr above the rms. The lower panels of 
Fig. 3 show the CO(l-O) wings amplitude binned in intervals of 
UV radius, covering baselines from 10 to 200 m. The decreasing 
visibility amplitudes are totally inconsistent with an unresolved 
source, which would instead give constant amplitudes with ra- 
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Fig. 2. CO(l-O) maps of the broad wings: the blue wing (left 
panel), integrated between -500 and -700 km s"' (right), and of 
the CO(l-O) red wing, integrated between 500 and 800 km s"' 
(right). The size of the maps is 15 "by 15". Each contour level 
corresponds to 0.2 mJy/beam. The beam size is 2.74"x 3.4". 
The cross indicates the position of the AGN. 




dius. The red symbols show the circular circular gaussian model 
fitted to the amplitudes. The fit is not perfect, suggesting that the 
real spatial distribution of the red and blue CO wings is more 
complex than a simple circular gaussian. We find that the red 
wing is spatially resolved with a significance of > 5cr while the 
blue wing is resolved only at 1.4cr. The FWHM of the com- 
bined red plus blue wings is 1.42"+0.2", corresponding to 1.2 
kpc (Table 1). The inclined disk model gives major and minor 
axes of 3.2" + 0.48"and 1.65" + 0.4", respectively, and a posi- 
tion angle of 68 + 1 1 degs. We note that the extended blue-shifted 
emission is consistent with the result on Na I D by Rupke et al. 
(2005). 

The determination of the outflow mass rate depends on the 
wind geometry and on the conversion factor, a, between the CO 
luminosity and the molecular gas mass. We estimated the CO lu- 
minosity by fitting the observed line profile with a narrow plus a 
broad component (see Fig. 1). The integrated CO luminosity of 
the broad component is L(CO)b - 1.16 x 10"^ K km/s pc^, about 
1/10 the luminosity of the narrow component. We converted the 
CO luminosity of the broad component into molecular gas mass 
M(H2) by assuming a conservative conversion factor a - 0.5 
M0 (K km s"' pc^)"', i.e. 1/10 the Galactic value . This is the 
lowest conversion factor found in different locations of M 82 (a 
typical starburst galaxy), including its molecular outflow (Weiss 
et al. 2001). We derive a mass of the outflowing molecular gas 
M(H2) - 5.8 X 10** Mq, which is consistent with the lower limit 
of 7 X 10^ Mq infen-ed by Fischer et al. (2010) based on the 
absorption molecular lines detected by Herschel. By assuming 
that this gas is uniformly distributed in a spherical volume of 0.6 
kpc in radius, its inferred density is ~ 25 cm"-*. Since the out- 
flow velocity is at least 700 km s"', the inferred mass outflow 
rate is dM{H2)/dt = 2200 Mq yr"'. If we assume density profile 
scaling as r^^ for the gas distribution, the inferred outflow rate 
is 710 Mq yr"' . An alternative, even more conservative estimate 
of the outflow rate can be derived by ignoring the Gaussian fit 
of the broad CO component and by using only the luminosity 
directly measured from the broad wings, i.e. by integrating their 
flux at velocities higher than -1-400 km s"' and lower than -400 
km s The luminosity of the wings is L(CO)= 3.2 x 10^ K km 
s"' pc^ (~30% of the CO luminosity of the total broad compo- 
nent inferred through the gaussian fitting). In this case we obtain 
a lower limit on the outflow rate of 600 M© yr ' , for the uniform 
gas distribution. By assuming a oc r^^ gas density profile we ob- 



Fig. 3. Residual maps of a point source fit to the blue (top left 
panel) and red (top right panel) CO wings. The cross indicates 
the position of the AGN. Amplitude (in Jy) versus baseline ra- 
dius for the blue (bottom left panel) and the red wings (bottom 
right panel). The red symbols show the results of a circular gaus- 
sian fit (see Tab.l). We recall that for an unresolved source the 
amplitude would be constant at all baselines. 

tain the most conservative lower limit is ~260 Mq yr"'. We note 
that alternative geometries of the outflowing wind, such as shell- 
like or disk-like configurations, give higher outflow mass rate, 
both because the inferred outflowing gas density is higher and 
because the de-projected outflow velocity is higher (e.g. if the 
molecular outflow occurs on the galaxy disk plane). Similarly, 
it is easy to show that in the case of a bipolar outflow the in- 
ferred outflow rate is the same or may be even higher (the line 
of sight has to intercept the outflowing molecular gas, since it is 
seen in absorption with Herschel, implying that the true, depro- 
jected radius of the bipolar outflow may be equal or larger than 
the projected size). In any case the inferred outflow rate is much 
larger than the star-formation rate in the host galaxy of 200 Mq 
yr"' (Taylor et al. 1999, Davies et al. 2004). The mass loss rate 
being much larger than the rate at which gas is converted into 
stars implies a phase of rapid quenching of star formation in the 
regions reached by the outflow (~1 kpc scale). 

4. Discussion and Conclusions 

The total amount of molecular gas in the galaxy disk as in- 
ferred from the integrated emission of the narrow component, 
by using the CO-to-H2 conversion factor a appropriate for 
ULIRGs (Solomon & Vanden Bout 2005), is M(//2) ~ 10'° 





velocity range 


L co(i-o) 


FWHM CO 




(km/s) 


K km/s pc^ 


(") 


Blue wing 


-500 -i- -700 


9.3e+07 


0.85(+0.61) 


Red wing 


500 H- 800 


2.3e+08 


1.56(+0.24) 


Red+Blue w. 


combined 


3.2e+08 


1.42(±0.2) 



Table 1. Integration limits, luminosity and uv fit results of CO 
blue and red wings. 
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Mq. Assuming the measured outflow rate of 710 yr the 
total molecular gas mass will be expelled from the host galaxy 
within 1.4x10^ years, implying the suppression of any signif- 
icant star-formation activity on the same timescale. Even with 
the more conservative outflow rate of 260 M© yr"' the timescale 
required to expel the molecular gas in the host galaxy would still 
be 4x10^ years. Such a timescale is shorter than the Salpeter 
time (~ 5 X 10^ years). We note that the quenching of star- 
formation has probably already started in the central regions. 
Indeed, stellar populations younger than 5 - 20 x 10^ years are 
not observed within the central kpc (Lipari et al. 2009). The 
total kinetic power of the out-flowing gas is 1.2x10'*"' ergs s"' 
(assuming an outflow rate of 700 km s"'), corresponding to a 
few percent of the AGN bolometric luminosity, LboI ~ 5 x 10'*^ 
ergs s"' (Lonsdale et al. 2003) which, for a black hole mass of 
~ 6 X 10^ Mq (Tacconi et al. 2002), corresponds to 6% of the 
Eddington luminosity. Such value of the kinetic energy is very 
close to that expected in the case of a shock wave produced by ra- 
diation pressure onto the interstellar medium (Lapi et al. 2005). 
The corresponding energy injected into the interstellar medium 
is Aii ~ fxLBoiXtACN ~ 5.5x10^^ erg, at least 4 orders of magni- 
tude larger than the overall contribution of Supernovae inferred 
from the observed stellar mass and age of Mrk 23 1 (Davies et 
al. 2004), and even larger than the energy injected by the ra- 
dio jets (Lonsdale et al. 2003). This strongly indicates the ra- 
diation field of the QSO as the primary engine at the origin of 
the observed outflow. Indeed, the observed outflow velocity of 
~750 km s"' implies a Mach number M ~ 30 - 50, indicat- 
ing a highly supersonic motion and, therefore, the formation of 
a shock front expanding in the interstellar medium of the host 
galaxy. Interestingly, this Mach number is in agreement with the 
expectation for the coupling between the QSO radiation field and 
the gas (Lapi et al. 2005). Indeed, the Mach number can be calcu- 
lated as ~ AE/E + 1, where E= KxTx Mgas/mp = 3 x 10^"* 
erg is the thermal energy of the interstellar gas (we assumed a 
temperature T= 4000 K and an interstellar gas mass Mgas ~ 10'" 
Mq within the shock), and AE = 5.5 x 10^^ erg is the over- 
all QSO energy dumped into the interstellar medium during the 
time taken by the shock to reach the observed radius of 0.6 kpc. 
The above computation yields M ~ 50, in good agreement with 
our estimate based on the measured outflow speed. This suggests 
a solution to the long standing problem of the transport of energy 
from the nucleus to the bulk of the galaxy: a highly supersonic 
shock transports outwards the energy accumulated in the centre 
by radiatively accelerated nuclear winds. Extended emission and 
kinematic signatures in the CO line profiles similar to those ob- 
served here are predicted by hydrodynamic simulations studying 
the effect of AGN feedback on the molecular gas (Narayanan et 
al. 2006, 2008). 

Due to the simultaneous presence of a strong wind, heavy 
X-ray absorption and of high star-formation rate, Mrk 23 1 has 
been regarded as one of the promising candidates of a QSO tran- 
siting from the obscured accretion phase, accompanied by vigor- 
ous star formation, where AGN feedback onto the host galaxy is 
in action, to the un-obscured phase (Page et al. 2004, Stevens et 
al. 2005) in the framework of the AGN-galaxy co-evolutionary 
sequence. These new observations of Mrk 231 provide one of 
the first direct observational evidences of QSO feedback that is 
dramatically affecting the evolution of its host galaxy. The QSO- 
driven giant molecular outflow is expected to expel the disk on 
short time scale, therefore halting star formation. This discovery, 
and in particular the inferred fraction of kinetic energy injected 
into the ISM relative to the QSO luminosity, confirms the expec- 
tations of models predicting a tight connection between the evo- 



lution of massive galaxies and the energy released by the accret- 
ing black holes harbored in their nuclei. The greatly improved 
sensitivity and expanded bandwidth of current and forthcoming 
millimeter interferometers will allow this kind of studies to be 
extended to larger samples of QSOs. Therefore, it will be pos- 
sible to directly verify whether the QSO feedback onto the host 
galaxy is really ubiquitous or not. 
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